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Abstract 

We consider small-scale spheroidal clusters of weakly 
interacting massive particles in our Galaxy as non- 
compact gravitational microlenses and predict the ap- 
pearance of caustics in the plane of a lensed source. The 
crossing of these caustics by a lensed star can produce a 
large variety of light curves, including some observed in 
actual microlensing events that have been interpreted 
as manifestations of binary gravitational lenses. We 
consider also observable effects during the gravitational 
microlensing of stars of non-zero angular size with a 
given brightness distribution across their disks by such 
an exotic objects as natural wormholes and objects 
whose space-time environment is described with the 
NUT metric. We demonstrate that, under certain con- 
ditions, the microlensing light curves, chromatic and 
polarizational effects due to the properties of the lens 
and the star disk brightness distributions can differ con- 
siderably from those observed for a Schwarzschild grav- 
itational lens, so that their analysis can facilitate the 
identification of such objects. 



in observations of gravitational microlensing of stars by 
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Various international teams - MACHO, 
PLANET, and others - have detected several 
hundred microlensing events for stars in the Magellanic 
Clouds and in the bulge of the Milky Way. However, 
the nature of the objects acting as gravitational lenses 
remains unclear. According t o the statistics of o bser- 



vations by the MACHO team (JAlcock et al.ll2000l ). the 
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1 Introduction 

The idea of detecting dark matter via its gravitational 
act ion on the light arri ving from dis t ant st ars suggested 
by iBvalko I (jl969l ) and iPaczvnski I (jl986l ) was realized 
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lens masses are in the stellar-mass range. In two cases, 
the masses exceed the limit for a neutron star, while the 
luminosi ties of these lenses a re no higher than that of 
the Sun (JBennett et al.ll2002l ): these objects are proba- 
bly black holes. 

A considerable fraction of the lenses may be dwarf 
stars i n our Galaxy. This is favored by statistic al argu- 
ments ( Romberg. Kompaneets. fc Lukash II 19951: Kerins 



■ 19971 ). and also by the analysis of light curves fB og- 
danov & Cherepashchuk Il998l : ICaudi et al. II2OO2I ). In 
two cases red dwarfs, that created the effect of gravita- 
tional lens, were directly observed by the Hubble Spac e 



dark matter- Telescope (lAlcock et al.ll200ll : lKozlowski et al. 1120071 ) 



As noted in the recent review of Kochanek. Schneider, 
& Wambsganss ( 20061 ) , the fraction of lenses from the 
dark matter can be sufficiently small. However, un- 
certainties in the properties of the gravitational lenses 
still admit a large variety objects of stellar masses as 
candidates for this role. 

Along with the black holes mentioned above, the 

possibilities of the lenses being small-scale clusters of 

weakly interacting massive particles (WIMPs) f Cure- 

yich. Zvbin. fc Sirota 19971: Sazhin. Yag ola. & Yakubov 

19961 : IZakharovlll999l : lBogdanovll200ll ) have also been 



considered. Lenses of more exotic types of matter have 
also been sug gested, in particular, tunnels in space-time 
(wormholes) ( Eiroa. Romero, fc Torres ll200lLSafonova, 
Torre s, fc Romero I2OO2I : iBogdanov fc Cherepashchuk 
I2OO2), as is predicted by general relativity. The papers 



Rahvar fc Nouri-Zonoz ( 20031) and Rahvar & Habibi 
(|2004[ ) considered the gravitational microlensing of 
point sources by hypothetical objects whose local re- 
gions of space-time are described with the NUT metric, 
and compared theoretical light curves to observations. 
The aim of the current paper is to study effects ob- 
served for gravitational microlensing of stars by the ex- 
otic objects. It is shown, in particular, that observa- 
tions of stellar microlensing evens can reveal fine pho- 
tometric, chromatic, and polarizational effects due to 
properties of the lens and lensed star: its angular size, 
the brightness and polarization distributions across the 
stellar disk, and the wavelength dependences of these 
characteristics. Though the manifestations of such ef- 
fects are fairly small, their analysis can contribute to 
identification of the lens type. 



2 Microlensing by Clusters of Weakly 
Interacting Massive Particles 

One approach to the solution of the dark-matter prob- 
lem supposes the existence of weakly interacting mas- 
sive particles (WIMPs) that formed at an early stage of 
the evolution of the Universe. Because of their gravita- 
tional interaction, these particles should have produced 
a hierarchical structure of clusters with certain scaling 
properties. The masses of the smallest objects are of 



the order of a solar mass. iGurevich. Zvbin. fc Sirota 



(|l997l ) fitted the mass density of the WIMPs in a small- 
scale object with the following function of the distance 
r from the center of the cluster: 



Po, < r < Tc 

p{r) = ^ pn{r/rc)", Vc < r < R^ 
0, r>R, 



(1) 



The exponent in (1) is also a scaling parameter, equal 
to a = 1.8. The radius of the core, which has a constant 
density po, is denoted r^ and can range from 0.05 to 
0.10 of the clus ter radius R^. According to Gurcvich. 



Zybin, & Sirota (J199l1 ). a typical cluster has a mass of 
M^ = 0.1 - 1.0 Mo and a radius of R^ ~ lO^"* - lO^^ 
cm. 

The particles interact weakly both with each other 
and with external matter, so that the detection of such 
clu sters seems verv unlikely. However. Gurevich & Zy- 
bin (|l995l ) pointed out that, if such objects are located 
along the line of sight connecting an observer and a 
normal star, they can deflect the light rays from the 
star; in view of the masses of the clusters involved, this 
can be considered gr avitational microlensing. Gure- 
vich, Zybin, & Sirota ( 19971 ) computed flux- variation 
curves for the gravitational microlensing of stars by 



non-compact lenses made up of the WIMP clusters. 
The shapes of these curves can be rather similar to 
those of a Schwarzschild lens, with a symmetric pro- 
file and a sharp maximum, albeit with broader wings. 
In their comparison of these curves with the first ob- 
servational data obt ained bv the OGLE group. Sazhin. 
Yagola, & Yakubov (|l996l ) showed that in two of eight 
cases the observed light curves can be better described 
by a non-compact lens model. Later, however, it was 
shown that the broadening of the wings of these light 
curves ca n also be explained in a simple r way - by either 
blending (jWozniak fc Paczvnski II1997T ) or by the radi- 
ation of the lens it self, if it is a normal star f Bogdanov 
fc Cherepash c huk Il998l ). 

Zakharov (jl999l ) analyzed a spherically symmetric 
nonsingular model for a non-compact gravitational mi- 
crolcns, which describes the expected density distribu- 
tion in a WIMP cluster better than earlier models. Un- 
der certain conditions, such a microlcns can have cir- 
cular caustic curve in the plane of a lensed source and 
produce a symmetric light curve with two maxima and 
a dip between them. Similar curves can also be pro- 
duced by binary gravitational lenses and have already 
been observed in some microlensing events. 

In reality, the spherical symmetry of the density dis- 
tribution in a non-compact object can break down. One 
factor that can lead to such a disruption of the symme- 
try is the tidal deformation of WIMP clusters as a re- 
sult of their gravitational interactions with each other 
or with stars of the Galaxy. The possibility of such 
defor mations was pointed out bv Gurevich. Zybin, & 
Sirota ( 19971 ). However, the probability of a close en- 
counter involving appreciable tidal forces is very low. 
Another factor - the oblateness of the particle cluster 
by rotation - seems more likely. 

Rotation is a general property of systems that form 
as a result of the gravitational contraction of matter. 
Rotation explains, for example, the observed oblate 
shapes of elliptical galaxies and globular star clusters. 
In the origin of rotation in such systems, an important 
role should be played by gas-dynamical effects, which 
are obviously absent in the case of weakly interacting 
particles. There is, however, a mechanism capable of 
inducing rotation of the particles exclusively via grav- 
itational interactions. This mechanism, first suggested 



bvlHovle I (Il949f ) and then developed in more detail by 
Peebles I ( 19691 ). is based on the mutual attraction of 
tidally deformed clusters. In this case, the orbital mo- 
mentum associated with the relative motion of the clus- 
ters is transformed into rotational momentum of each 
of the interacting clusters. In the early Universe, when 
WIMP clusters formed, their spatial density was sub- 
stantially higher than now. Close encounters occurred 
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rather frequently, and the Hoyle - Peebles mechanism 
could be fahdy efficient. 

The observed manifestations of a transparent gravi- 
tational lens depend on the distribution of the surface 
mass density S(^), where ^ is a vector defining the posi- 
tion of a point in the plane of the lens. Here and below, 
we write vector names in bold font. For the spherically 
symmetric law (1) in a coordinate system whose origin 
coincides with the projection of the center of the parti- 
cle cluster onto the plane of the lens, S(f) is given by 
the integral 



^/R^^ 



S(0 



p{^e + ^) dz. 



(2) 



It is obvious that the distribution of the surface mass 
density in this case is rotationally symmetric. Further, 
we consider the dimensionless density 



k{0 = S(C)/Sc 



(3) 



where the so-called critical density Ecr is determined 
by the relation 



c'D. 



AnGDdDds 



(4) 



In (4), c is the speed of light, G the gravitational con- 
stant, Ds the distance between the observer and lenscd 
source, Dd the distance between the observer and lens, 
and Dds the distance between the lens and source. 

Let us adopt the radius of the cluster R^ as a scaling 
parameter ^o and introduce dimensionless vectors in the 
plane of the lens x = ^/^o and the plane of the source 
y = 77/7/0, where r/o = ^oDg/Dd, Co = Rx, x(a;i,X2), 
and 7(7/1,2/2)- We use a(x) to denote the normalized 
vector a*'(Cox): 



a(x) = 



DM 



ds 



^oDs 



aO(Cox). 



The gravitation al lens equation fKochan ek, Schneider, 
& Wambsganss uOOa : IZakharov I Il997l ) can then be 
written in the form 



y = x-a(x). 



(5) 



The possibility of WIMP clusters being oblate by 
rotation can hardly provoke any fundamental objec- 
tions. However, the questions of the oblate cluster's 
shape and density distribution have no simple answers. 
The branch of celestial mechanics that studies the equi- 
librium shapes of rotating bodies assumes them to be 
comprised of a viscous incompressible fluid, which is 
clearly inappropriate for systems of weakly interacting 



particles. Nonetheless, some conclusions may be of in- 
terest for our case as well. Given the presence of a 
small dense core and a rapid decrease of the particle 
density (1) with distance from the cluster center, we 
can describe this density distribution to a first approxi- 
mation using the Roche model, which assumes that the 
entire mass of the rotat ing body is concentrated at its 
center (jSubbotin Ill949l ). This type of model has been 



used earlier to study planetary atmospheres and proto- 
planetary disks. The oblateness of the marginally sta- 
ble surface of the Roche model was fairly small, with 
a polar-to-equatorial-axis ratio equal to 2/3 f Subbotin 
' 194911. 



Models with a nonsingular power-law elliptical dis- 
tribution of the surface mass density have been widely 
used in studies of transparent gravitational lenses: 



k(x) 



q 



(7/2 



4/k'y 



(6) 



This model has a constant-density core of radius u at 
the center. The contours of equal surface density are el- 
lipses with fixed minor to major semiaxis ratio k = b/a. 
Models of the form (6) are believed to adequately de- 
scribe the surface mass density distributions in galaxies 
for a wide range of degrees of oblateness. 

We first analyzed the case of a non-rotating cluster 
and tried to fit the density distribution of a WIMP clus- 
ter (3) with Vc = O.OlRx and normalization k(0) = 1 us- 
ing the model (6) with /s = 1, choosing model parame- 
ters to minimize the sum of squared residuals. The solu- 
tion of the corresponding optimization problem yielded 
the best-fit parameter values q = 1.04, u = 0.09, and 
n = 0.68. The surface density distribution for model (6) 
with these parameters is a good fit to the initial distri- 
bution. We then assumed that, as the particle cluster 
rotates, the parameters q, u, and n remain constant 
and the ellipses of equal density experience oblateness 
equal to that of the marginally stable surface in the 
Roche model: k = 2/3. 

It is well known that a transparent gravitational lens 
with surface density distribution (6) can have caustics 
in the source plane, which, when crossed by a lensed 
star, can produce a large variety of light curve shapes. 
For this to happen, it is sufficient for a point x with 
k(x) > 1 to exist in the plane of the lens fKoc hanek. 
Schneider, fc Wambsganss l2006l: IZakharovlll997h . This 
condition can be satisfied by the WIMP clusters consid- 
ered. Let us fix Rx = 10" cm and A'4 = 0.5 Mq = 10^^ 
g. It then follows from (1) that po = 1-8 x XQ'^gjcnv' , 
and the computation of integral (2) yields a surface 
density of I](0) = 3.4 x Wf'glcrr? at the lens cen- 
ter. We find from (4) that the critical density is 
Scr = 5.3 X lOVcTO^ and Scr = 1-7 x lO^g/cm^ for 




Fig. 1 Caustic curves of tlie non-compact WIMP lens in 
tiie plane of the source 



typical cases of Icnsing of stars in the Large Magellanic 
Cloud by an object located in the Galactic halo {Dd = 8 
kpc, Dg = 46 kpc, and Dds = 38 kpc) and lensing of 
Galactic bulge stars {Dd = 4 kpc, Dg = 8 kpc, and 
Dds = 4 kpc), respectively. Thus, the required excess 
of the density over its critical level is quite possible. 

Let us assume for the sake of simplicity that the ro- 
tational axis of the cluster lies in the plane of the lens 
and that the major axes of the ellipses of equal surface 
density form intervals on the xi axis. It is necessary 
however to note that this alignment of cluster maxi- 
mizes th e asymme t ry of the lens and the size of the 
caustics. 



Barkana I (|l998l ) has developed an efficient al- 



gorithm for solving the gravitational-lens equation (5) 
using the density distribution (6) and written a FOR- 
TRAN code implementing this algorithm. For a given 
vector X, the code computes the components of the de- 
flection angle vector a(x) and the elements of the Ja- 
cobi matrix Aij = dyi/dxj . Our computations for 
model (6) with the parameters adopted here demon- 
strated the existence of critical curves in the plane of 
the lens, defined by the condition detA = 0. The corre- 
sponding caustics in the plane of the source are shown 
by the solid curves in Figure 1. When a point source 
crosses the caustic, the gravitational lens amplification 
factor fi = l/|(ietA| can formally approach infinity. For 
a source with a finite angular size, /i is finite but can 
achieve rather large values. 

For the adopted size of the WIMP cluster, the unit 
scale length in the plane of the lens for Dd = 4 kpc cor- 
responds to an angle of 0".0017, which is much smaller 
than the resolution of current telescopes. As a result. 



the object should show up as a microlens, with the total 
flux from all its images being observed. Due to the or- 
bital motion of the observer, particle cluster, and lensed 
star, the observed flux should vary with time. Exact 
computation of the light curves is computer intensive. 
However, it is easy to obtain a qualitative description of 
the observational picture. When a source of finite angu- 
lar size crosses the outer caustic, two additional source 
images should appear and the observed flux should in- 
crease. However, the flux amplification factor on this 
caustic is comparatively small, and the caustic cross- 
ing could escape detection. When the source crosses 
the inner caustic, two more source images appear, re- 
sulting in a sharp increase of the flux followed by a 
slower decrease. At the second crossing of the inner 
caustic, the two images merge and disappear, inverting 
the pattern of the flux variations. When the line of 
the relative motion of the lensed star in Fig. 1 passes 
far from the cusp and is parallel to the coordinate axes 
when it crosses the inner caustic, the light curves have 
symmetric profiles with two maxima and a dip at the 
center. Such curves can be re alized in the nonc ompact 
microlens model analyzed by IZakharov I (jl999r ) . If the 



star moves along a line at some angle to the yi axis, 
either the first or second crossing of the inner caustic 
should take place closer to the corresponding cusp. As a 
result of the difference in the amplification factors, the 
light curve becomes asymmetric and either its first or 
second maximum becomes stronger. Similar light-curve 
shapes can appear in the cases of binary gravitational 
lenses and have been observed in actual microlensing 
events. We must therefore bear in mind the possibility 
that at least some such events could in reality be associ- 
ated with WIMP clusters. When the star goes beyond 
the outer caustic, comparative modest fiux variations 
will again be observed as a flux increase followed by a 
more abrupt decrease. If the lensed star does not cross 
the outer caustic, the flux-variation curve is symmetri- 
cal and has a single maximum. 

The microlensing pattern considered above depends 
on a large number of parameters. Obviously, it is pos- 
sible to obtain a large variety of fight curves by chang- 
ing the lens orientation and varying the parameters for 
the density distribution and the relative motion of the 
lensed star. Represents the interest to examine quanti- 
tative a question how well a light curve for binary lens 
can be fitted by curve for the model of noncompact lens. 



3 Spatial - Temporal Tunnels as Gravitational 
Microlenses 



As was s hown by Morris fc Thorne ( 19881). Hochberg 
& Visser ( 19971 ). and Hochberg, Popov, & Sushkov 
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([1993)1 the equations of general relativity allow solu- 
tions in the form of two regions in space-time connected 
by a so-called spatial-temporal tunnel or natural worm- 
hole. Such a tunnel can exist only if it is filled with 
a certain type of exotic matter with a negative energy 
density a nd, therefore, a nega tive mass. It was first sug- 
gested bv lKim fc Choi (|l996l ) that such an object could 
be manifest as a gravitational lens capable of forming 
multiple images of a distant source of radiation. Under 
certain conditions, the angular distance between the 
images can become smaller than the resolving power of 
telescopes, and the observer will detect only the com- 
bined flux from all the images of the source. In this 
way, the wormhole will be manifest as a gravitational 
microlens with negative mass. Such a microlens mod- 
eled as a negative p oint mass was considered bv Eiroa. 
Romero, fc Torres (l200ll'). and Safonova. Torres, & 
Romero (J2002l ) using the commonly applied approxi- 
mations of geometric optics and small deflection angles. 
Like an ordinary gravitational lens, it is basically speci- 
fied by the angular radius of the Einstein cone poi which 
can be expressed as 



pI 



AG\M\ 



D, 



{Dds + Dd)Dd 



(7) 



where M the mass of the lens. Unfortunately, any re- 
liable estimations of mass of wormholes and hence po 
until today are absent. 

The lens equation indicates that its properties differ 
dramatically from those of an ordinary positive- mass 
Schwarzschild lens due to the presence of a circular 
caustic with angular radius 2po in the source plane. 
When a point source is located an angular distance 
p < 2pQ from the lens, it is no longer observed, and 
the detected flux decreases to zero. When p > 2po, 
two images of the source arise, located on one side of 
the lens between the directions towards the lens and 
the source. One of the images is always inside the Ein- 
stein cone, while the second is outside it. As the source 
approaches the caustic from the outer side, the images 
merge and disappear when p = 2pQ. The ampliflcation 
factor for the flux from a source with an inflnitely small 
angular size in the region p > 2po is 



A{u)^ 



lyi? 



(8) 



where il = p/pa] at a caustic crossing, A{u) formally 
tends to infinity. Relation (8) differs from the corre- 
sponding formula for a Sc hwarzschild lens fKoch anek. 
Schneider, & Wambsganss 



warzschild lens fKoch anek , 
200dlZakharovlll997l ) only 



in the signs of the numerical terms in the numerator 
and denominator. 



In the case of lensing of a source with a finite angular 
size, the flux will increase as it approaches the caustic, 
then sharply decrease and completely disappear. In the 
course of the egress of the source from the inner region 
of the circular caustic, the pattern will change to the 
opposite. For a source with a circularly symmetrical 
brightness distribution, the light curve will always be 
symmetrical relative to the time of its closest approach 
to the lens. If the source does not cross the caustic 
during its relative motion, the light curve displays one 
maximum and its shape is close to those of observed 
stellar microlensing light curves. A s it was shown by 
Bogdanov fc Cherepashchuk I (J2002r ). if a point source 
moves outside the circular caustic, photometric obser- 
vations cannot distinguish light curves due to a worm- 
hole and a positive-mass gravitational lens. However, it 
is apparent that the ampliflcation factors of the lenses 
can be equal only for one value of the impact param- 
eter, so that the light curves for extended sources will 
differ. 

Images and light curves of extended sources formed 
by a negative-mass lens were calculated bv Safonova. 
Torres, fc Romero ( 2002f ) using the method of ray 
tracing and the solution of the lens equation on a 
5000 X 5000 count grid. If we consider only photomet- 
ric effects, the problem is simplifled by the fact that 
the ampliflcation factor for the flux from a source area 
element will be specified by (8). The problem is thus 
similar to the calculation of light curves for the case 
when a source is l ensed bv a linear caustic ( Bogdanov 
fc Cherepashchuk |2000|). In this case, photometric ef- 
fects depend only on the one-dimensional projection of 
the source brightness distribution onto the axis perpen- 
dicular to the caustic usually called the strip brightness 
distribution. In our case, the lens possesses a circular 
caustic, and the light curve I{p) as a function of the an- 
gular distance from the lens p will depend on the one- 
dimensional brightness distribution Bp{r), which is the 
source brightness integrated over an inflnitely narrow 
circular strip with radius r concentric to the caustic. 
Within this strip, the lens amplification factor (8) re- 
mains constant. If the angular size of the source is much 
smaller than pq, the curvature of the caustic can be ne- 
glected, and the distribution Bp(r) coincides with the 
strip brightness distribution. As in the previous case 
of a linear caustic, the influence of the singularity of 
the integral in the flux calculation can be removed by 
corresp onding selection of the grids for the variables pi 
and ri (JBogdanov fc Cherepashchuk II2002I ). 

Further, we will assume the lensed source to be a star 
and the brightness distribution over the stellar disk to 
display circular symmetry. In this case, the detected 
flux will depend only on the angular distance between 



the center of the stellar disk and the lens. For most 
stars, the limb darkening is described by the linear law: 



6a (a*) = bx{l){l -x\+ x\y), 



(9) 



where /i is the cosine of the angle between the line of 
sight and the normal to the stellar surface and x\ is 
the limb-darkening coefficient, which depends on the 
wavelength A. The value ^ corresponding to the angular 
distance from the center of the star p is related to it by 
the expression; 



/i = V^R23^/p^ 

where R is the angular radius of the star. Further, we 
will assume that the brightness at the center of the star 
can be written 

6a(1) = lAi?'(l - xa/S), 

which normalizes its flux to unity. 

Let us introduce a polar coordinate system (r, if) 
with its origin coincident with the lens. If the center 
of the stellar disk is located an angular distance p from 
the lens, the integrated brightness of the circular strip 
of the disk specified by the radius r is 



Bp{r) 



bx{p)r dip, 



(10) 



where p^ ^ r^ -\- p^ — 2rp cos (p and an upper limit for 
the integral can be derived from the condition R^ = 
^2 _|_ p2 _ 2^p QQg ip^^ Taking into account the caustic 
amplification factor (8), the detected radiation flux is 



P+B. 



I{P) 



A{r/po)Bp{r) dr, 



(11) 



p-R 



where A{r/pQ) = for r < 2po . 

For the purpose of the estimation of the effect of 
finite angular size of the source we chose a KOV star, 
and assumed that its l imb-darkening coeffic ient is xy = 
0.702 in the V band |Rubashevskii Ill990l ). The light 
curves were calculated using (9), (10) and (11) and 
canonical-dissection grids with a step sufficiently small 
that the relative error would not exceed 10""'. Since 
the flux depends only on the angular distance p, the 
light curve can be calculated for a single value of the 
impact parameter pm = 0; we can find the counts for 
other values via interpolation. We used a cubic spline 
interpolation, which provided the necessary accuracy. 
We chose the value R = 0.10 for the comparison with 



the calculations of lEiroa. Romero, fc Torres I (|200lD 




t/t. 



Fig. 2 V light curves for the microlensing of a star with 
angular radius 0.10 by a wormhole; time is measured in 
units of the Einstein cone angular radius crossing time. The 
positive branches of the curves for the impact parameters 
Pm = 2.1, 1.9, 1.0, and are indicated 

The light curves were calculated for the impact pa- 
rameters Pm ~ 0,1.0,1.9, and 2.1. Further, we will 
measure time in units of the time required to cross 
the angular radius of the Einstein cone tE ■ Figure 2 
presents these results. Since all the curves arc symmet- 
rical about the time of closest approach to the lens, the 
figure shows only their positive branches. Our calcula- 
tions are consistent with the results of Eiroa. Romero, 
& Torres ( 200ll ). confirming the validity of the tech- 
nique used. We can see from Fig. 2 that, when the 
star does not cross the circular caustic, the light curve 
displays a single maximum and its shape is similar to 
that for a Schwarzschild lens. 

When a star crosses the circular caustic, character- 
istic symmetrical light curves with a central depression 
arc formed. They differ markedly from observed stellar 
microlensing light curves, and could potentially provide 
unambiguous evidence for the presence of a negative- 
mass gravitational lens. However, in the studv bv Bozza 



(further, we will measure all angles in units of po)- 



et al. ( 2002) of microlensing of stars by a compact ob- 
ject surrounded by an extended gaseous envelope, due 
to the extinction of the stellar radiation by refraction 
and Rayleigh scattering in the envelope, the light curve 
can acquire a shape close to that characteristic of a 
negative-mass lens. One possibility for discriminating 
between these two types of lens is provided by analyses 
of chromatic effects associated with the lensing. 

The basic property of gravitational lenses of both 
positive and negative masses is that their parameters 
are independent of wavelength. However, the amplifi- 
cation factors are different for radiation from different 
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regions of the stellar disk. Due to the wavelength de- 
pendence of the brightness distribution, this results in 
chromatic effects during the microlensing. The possibil- 
ity of observing these effects for ordinary gravitational 
lenses has been considered, for exa n iple, in papers of 
Bogdanov fc Cherepashchuk I (119951'). iLoeb fc Sasselov 



(|l995l ). and lHan. Park, fc Jeongl |2000l ). Such an anal 
ys is for a negative-mass lens was c arried out in paper 
of lEiroa. Romero, fc Torres I (|200ll ) by calculating the 
U and / light curves for a K giant with T^ff = 4750 K 
and R = 0.10, assuming that the limb-darkening co- 
efRcients were xjj = 1.050 and xj = 0.503. The limb 
darkening coefficients x\ were determined from calcula- 
tions of stellar atmosphere models. In this case, xx can 
be estimated in different ways, based on either conser- 
vation of the total flux or th e best consistency be tween 



(9) and the calculated data (JRubashevskii 



Il990h . The 

first method can formally give a value a;A > 1 ■ How- 
ever, it is physically meaningless to use this value to 
calculate t he brightness distribution (9). a s was done 
in paper of lEiroa. Romero, fc Torres I (J200l( ). since this 
yields negative brightness at the edge of the stellar disk. 

We calculated the variations oiV — I for microlens- 
ing of a KOV star by a negative-mass object, since ex- 
actly these photometric bands are most frequently used 
in searches for gravitational- microlensing events. The 
limb-darkening coefficient s xy — 0.702 and xj ~ 0.433 
were taken from paper of iRubashevskii I ( 1990|). As for 
the calculation of the photometric effects, we adopted 
R = 0.10 andp™ = 0, 1.0, 1.9, 2.1. Figure 3 presents the 
resulting positive branches of the color- variation curves. 
As expected, the maximum variation A(V — /), which 
corresponds to reddening, occurs when the stellar disk 
is almost entirely inside the circular caustic. However, 
the radiation flux is very small in this case (see Fig. 2) 
and cannot be detected. We restricted the calculations 
to the region in which the V flux exceeds 0.001 of its ini- 
tial value. In the course of the passage of the stellar disk 
across the caustic, A{V — I) changes sign twice: after 
reddening, the star becomes bluer, then redder again. 
However, in the region where the flux is appreciably 
different from zero, the chromatic effects are small and 
can be detected only with rather high photometric ac- 
curacy. If the star does not cross the circular caustic, 
the c olor varies in the same way as for a Schwarzschild 
lens (JBogdanov fc Cherepashchuk Ill995l ) and reddening 
is always observed. 

Thus, in the same way as for ordinary gravitational 
lenses, chromatic effects during the microlensing of a 
star by a wormhole are second-order compared to pho- 
tometric effects. However, their detection could make 
it easy to discriminate between this exotic object and a 
compact body with an extended gaseous envelope. In 
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Fig. 3 Variation of the color index A{V — I) for tlie 
microlensing of a star by the wormhole shown in Fig. 2. 
The units along the horizontal axis are the same as in Fig. 
2 



the latter case, the chromatic effects will be very differ- 
ent in both their values and character, due to the strong 
wavelength dependence of the refraction and, especially, 
of Rayleigh scattering (oc A~*). 

Analogous to chromatic effects, gravitational mi- 
crolensing of stars gives rise to variable polarization of 
the observed radiation, due to the character of the po- 
larization distribution over the stellar disk. The stellar 
radiation becomes polariz ed as a result o f scatt ering 
in the stellar a tmosphere. IChandrasekhar I (jl953r ) and 
SobolevI ( 19561 ) demonstrated that the radiation emerg- 
ing from a plane-parallel Rayleigh atmosphere is par- 
tially polarized, with the direction of the electric vector 
being perpendicular to the plane containing the line of 
sight and the normal to the stellar surface. In ordinary 
observations of stars, the detected total radiation flux 
from the entire stellar disk is unpolarized due to the 
disk's symmetry. The origination of partial polariza- 
tion of the total flux is associated with a violation of 
the disk's symmetry due to rapid rotation of the star, 
tidal deformation, the presence of hot spots and/or an 
asymmetrical envelope, or a combination of these fac- 
tors. It is apparent that the action of a gravitational 
lens also violates the disk's symmetry. The type and 
degree of polarization observed during microlensing of 
a star by a Schwarzschild lens we r e cons idered in papers 
of lSimmons. Willis, fc Newsam I (1995) and Bogdanov. 
Cherepashchuk fc Sazhin (J19961 ). 

To analyze the linear polarization that originates 
during the propagation of radiation in the atmosphere 
of a star, it is convenient to introduce the intensities 
bi in the plane, which contains normal to the stellar 



surface and line of sight, and br in the direction per- 
pendicular to this plane. For simplicity, we will assume 
that the lensed star has a plane-parallel Rayleigh at- 
mosphere. In this case, the exact brightness distribu- 
tions bi and br o ver th e stellar disk are presented by 
(Il953h . 



4,0 r P{t) 



Chandrasekhar 



For our estimates, it is more 
convenient to use the app r oxima te formulas derived by 
Bochkarev fc Karitskava I (j 19831 ): 



hr + bi = 



1 + 16.035 ^i + 25.503 ^^ 
1 + 12.561 ^-1-0.331 ^2 



0.1171 + 3.3207^^ + 6.1522/^2 



(12) 



(13) 



l-^l 1 + 31.4160 /i + 74.0112 /i2 ' 

where, as in (9), /i is the cosine of the angle between the 
line of sight and the normal to the stellar surface. Re- 
lations (12) and (13) can be used to estimate the inten- 
sities of the polarized components, so that the relative 
error of the degree of polarization P = {br^bi)/{br + bi) 
compared to the exact solution for any /i does not ex- 
ceed 0.001. 

Let us bring the origin of the polar coordinate sys- 
tem (p, 6*) to the center of the stellar disk, so that the 
main axis from which we measure the angle 9 will be 
directed towards the gravitational lens, located an an- 
gular distance p from the coordinate origin. Due to 
the circular symmetry of the initial brightness distribu- 
tions bi and br, it is apparent that the flux and degree 
of polarization of the detected radiation will depend 
only on p. Let bL{p,6) and bfi{p,9) be the intensities 
at the point of the visible stellar disk with the coordi- 
nates (p, 9) for observations through polarizers oriented 
parallel and perpendicular to the direction towards the 
lens. Then, according to the transformation law for the 
Stokes parameters (and, accordingly, for the intensities 
bi and br), we have: 



6l(p, 9) = br{p) sin^ 9 + bi{p) cos^ 



bji{p, 9) = bi{p) sin^ 9 + br(p) cos^ 



(14) 



(15) 



The coordinate system [p, 9) is related to the system 
(r, (/)) introduced above: rsintp = psin9. Taking this 
into account and substituting the brightness distribu- 
tions (14) and (15) into (10) and (11), we can calcu- 
late the flux from the star with the polarizer oriented 
in the parallel (/^(p)) and perpendicular {Iji{p)) direc- 
tions relative to the lens. The degree of polarization 
P{p) when the center of the star is an angular distance 
p from the lens will be 



Pip) 



In{p)-h{p) 
Inip) + Il(.p)' 





(16) 



t/t. 



Fig. 4 The variation of the degree of polarization P{t) (in 
percent) for curves of microlensing of a star by the wormhole 
shown in Fig. 2. The units along the horizontal axis are the 
same as in Fig. 2 



Figure 4 presents the positive branches of the po- 
larization curves for the microlensing of a star with a 
Rayleigh atmosphere calculated for the same lensing pa- 
rameters as in the previous case. As previously, we will 
consider only those sections of the curves for which the 
flux is no lower than 0.001 of the initial value. Since, 
like the chromatic effects, the polarization effects of mi- 
crolensing are due to differences in the amplification of 
the radiation from the central and limb regions of the 
stellar disk, it is not surprising that the shape of the po- 
larization curves essentially coincides with that of the 
color-index variation curves (see Fig. 3). During the 
passage of the caustic across the stellar disk, the degree 
of polarization changes sign twice. According to (16), 
this means that the plane of the electric vector is first 
perpendicular {P > 0), then parallel (P < 0), and fi- 
nally perpendicular {P > 0) to the direction towards 
the lens. 

When the star does not cross the circular caustic, 
the degree of polarization is always positive, and the 
plane of the electric vector is always perpendicular to 
the direction from the center of the star towards the 
gravitational lens. Thus, with sufficient observational 
accuracy, it may be possible to detect the rotation of 
the plane of polarization as the star moves relative to 
the lens. At the time when the maximum degree of 
polarization is attained, which coincides with the max- 
imum of the flux from the lensed star, the plane of po- 
larization will be specified by the position angle of the 
relative motion of the star and lens. The direction of 
the motion remains unspecified, due to the symmetry. 
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A similar pattern is observed for a Schwarz schild lens 
i Bogdanov. Cherepashchuk fc Sazhin Ill996l ). 

For a star with a pure Rayleigh atmosphere, the de- 
gree of polarization is small, and its maximum, reached 
at the very edge of the stellar disk, is 11.7%. Since the 
albedo of a singular scattering under the conditions in 
stellar atmospheres is always smaller than unity, we ex- 
pect that the polarizations of real stars will be lower, so 
that the calculated values correspond to upper limits. 
In addition, the Nagirner effect may also be observed 



(JNagirner I Il962l ) , which is manifest by a rotation of 



the plane of polarization with motion from the edge of 
the stellar disk. The influence of the true absorption 
on the polarization of the radiation emerging from the 
atmosphere of the star assuming a plane-parallel geom- 
etry was consi dered in detail bv Bochkarcv. Karitskaya, 
& SakhibuUin ( 19851 ) . In the Wien region of the spec- 
trum, P can appreciably exceed the value for a Rayleigh 
atmosphere. The degree of polarization for long wave- 
lengths is smaller, but remains appreciable. For late- 
type dwarfs, which statistically are most likely to be 
lensed, an additional contribution to the increase of P 
may be introduced by Rayleigh scattering on molecules. 
Overall, for real stars, the order of magnitude of the 
calculated polarizations is expected to remain valid. 

Thus, both chromatic and polarization effects dur- 
ing the microlensing of a star by a wormhole are small. 
However, the accuracy of current polarimetric obser- 
vations considerably exceeds photometric accuracy; for 
sufficiently bright objects, P can be determined with 
an error not exceeding 0.01 %. Unfortunately, this 
accuracy cannot be provided in monitoring observa- 
tions during searches for microlensing events. How- 
ever, observers have already accumulated considerable 
experience in triggering large telescopes to observe de- 
tected events involving microlensing by gravitational- 
lens caustics, so that prospects for observing these po- 
larization effects are fairly encouraging. 



4 Microlensing of Stars by NUT 
Gravimagnetic Lenses 

Newman, Unti, and Tamburino (NUT, 1963) consid- 
ered a space-time metric that is a generalization of the 
Schwarzschild metric. Along with the mass of ordinary 
matter, Af , it is also determined by the so-called NUT 
factor, I, which describes the contribution of magnetic 
monopole to the curvature of space-time. The parame- 
ter I is often simply called the magnetic mass. The NUT 
metric becomes the Schwarzschild metric for I -^ 0. On 
the other hand, gravitation still exists as M -^ (the 
"pure NUT" case). 



The ability of the N UT objects to be gravitational 



lenses was studied by iNouri-Zonoz fc Lvnden-Bell 
( 19971 ) . Generally, the properties of such lenses are close 



to those of Schwarzschild lenses. Under certain condi- 
tions, they are also capable of forming multiple images 
of lensed objects and amplifying the light flux detected 
by the observer. However, there are two important dif- 
ferences: the presence of a round zone where the lensed 
source is invisible, with its center towards the lens and 
its size increasing with the NUT factor, /; and a shear of 
the source images, manifests as a rotation in the plane 
of the sky around the optical axis of the lens. 

It is possible that the angular separation between 
the source images formed by a NUT lens will be smaller 
than a telescope's resolution. In this case, we are deal- 
ing with gravitational microlensing, and the observer 
will detect only variations of the light flux during the 
lens's motion relative to the source. This situation was 
analyzed for NUT lensin g of a point so urce in papers of 
Rahvar &: Nouri-Zonoz I (|2003i ) and iRahvar fc Habibi 



(200J), where the results were also compared to lensing 
by a Schwarzschild lens. 

The amplification coefficient for the flux from a 

i ioint source lensed by a Schwarzschild lens i s given by 
Kochanek. Schneider, fc Wambseanss II2006I: Zakharov 
19971 '! 



A{u) 



u\/v? 



(17) 



Here as earlier, u ~ p/po is the angular distance be- 
tween the lens and the source measured in Einstein cone 
angular radii. When describing microlensing by a NUT 
object, a new parameter is introduced: 



Po V DdDs 



When u^ < 2{^/4:N'^ + 1 - 1) , the amphfication co- 
efficient of a NUT lens is imaginary, corresponding to 
an absence of source images. Outside this region, the 
amplification coefficient for the flux from an infinitely 
small source is given bv lRahvar fc Nouri-Zonoz I (J2003l ) 



A{u) 



8JV^ 



-1/2 



8W* 



-1/2 



(18) 



where 



a± 



±V^ 



4u2 - 167V4 



2(47V4 _^ 1) 



and the sign in the subscript of a corresponds to the sign 
in front of the square root. An approximate expression 
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for the amplification coefficient valid for small N and 
A{u) is 



A{u) 



u2 + 2 SN'^iu^ + 2) 

iV^?T4'^w3(u2+ 4)3/2- 



(19) 



Comparing (f9) and (f7), we note that amplification 
by a NUT lens is stronger than amplification by a 
Schwarzschild lens when m^ > 2{^/AN^ + 1- I) . With 
decreasing magnetic mass (N — > 0), the second term in 
(19) vanishes, and this equation becomes equal to (17). 
The presence of the source invisibility zone for N not 
equal to zero limits the minimum value of u, and de- 
creases the maximum possible amplification coefficient 
for a NUT lens, compared to a Schwarzschild lens. 

In a microlensing event for a real star with finite an- 
gular size, the flux from each infinitesimal area of its 
disk will be multiplied by the amplification coefficient 
(17) or (19), depending on the lens type. Let us assume 
for simplicity that the brightness distribution across the 
stellar disk, b\, which depends on the wavelength A, can 
be represented with a linear law (9). Let r be the stel- 
lar angular radius and p the angular distance from the 
stellar disk center measured in Einstein angular radii 
Pq. The brightness distribution (9) can then be written 



b\{p) ^bl{l^xx+ xxy'r^-p'^/p), 



(20) 



where &° the brightness at the center of the stellar disk. 
Searches for microlensing events are based on stellar 
photometry in a system close to the standard broad- 
band UBVRI system. We can write expressions similar 
to (9) and (20) for the brightness distributions across a 
stellar disk in these filters, introducing the brightness at 
the disk center and the limb-darkening coefficient for a 
given filter, for example, by and xy for the V filter. We 
obtain the flux in this fllter in the absence of the lens, 
Hy , by integrating the brightness distribution over the 
entire visible stellar disk: 



Hy^Trr^b1.{l-xy/3). 



(21) 



The limb-darkening coefficients in the filters of the 
standar d broadband svstem were computed bv Ruba- 
shevskii ( 1990) for stars with various effective temper- 
atures using Kurucz' grid of model atmospheres. 

In the presence of a gravitational lens at an angular 
distance u from the center of the stellar disk, the V flux 
can be written 

27r r 



Hyiu) = J d^ J A{^u^ + e ' ^u^ COS ^)x 

feO.(l -XV + xv^R^^eir)i d^, 

where the angle ip is measured from the direction to- 
wards the lens and the ampliflcation coefRcient A is 



given by (17) for a Schwarzschild lens or (18) for a 
NUT lens. Similar expressions can also be written for 
the fluxes in other filters. Taking into account the de- 
pendence of u on time, t, 




where uq is the impact parameter - the minimum angu- 
lar distance, achieved at the moment to, and ts is the 
time needed to cross the Einstein cone angular radius, 
we find the observed light curve, /y(i): 



Iy{t)=Hy{t)/H^. 



(23) 



As where mentioned in section 3, the amplification 
coefficient for any gravitational lens is independent of 
wavelength. However, in the most probable case when 
the lens is not projected onto the stellar disk, the rapid 
decrease of A{u) with increasing distance, u, leads to 
an increase in the fraction of light from the disk limb 
in the observed flux. For the UBVRI filters, the limb 
darkening decreases with increasing maximum sensitiv- 
ity wavelength. Thus, a chromatic effect, a kind of red- 
dening, should be observed for real stellar microlensing 
events. This chromatic microlensing effect can be re- 
vealed by measuring color indices, for example, 

B-V^2.5 \g{Hy/HB) + Cb-v 



V^R^2.5lg{HR/Hy 



C 



V-R, 



with the constants Cb-v and Cy-R determined by the 
zero points of the photometric system. In the case of 
reddened star light, the value of the color indices will 
exceed the values in the absence of the lens, {B — V)o 
and {V — R)o, and should be largest when the lensed 
star achieves its maximum light flux. 

We computed the photometric and chromatic effects 
of microl ensing a star by a NUT l ens as suming, as in 
paper of iRahvar fc Nouri-Zonoz I ( 20031) and Rahvar 
& Habibi ([20041), that < iV < 0.5. We chose a 
main-sequence KOV star with an effective temperature 
of Teff ^ 4900 K and color ind ices {B - V)o = 0.890, 



{V - R) n = 740 (lAllenlll973l). According to Ruba- 



shevskii ( 1990) , the limb-darkening coefficients for such 
a star derived from the condition of flux conservation 
are xb = 0.859, xy = 0.702, and xr = 0.585. We se- 
lected the brightnesses at the disk center so that the 
computed fluxes (21) in the absence of the lens would 
give the original color indices, {B — V)o and {V — R)o. 
We assumed the most probable case of lensing, when 
the NUT- lens invisibility zone does not overlap with 
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Fig. 5 Light curves for microlensing of a star with angular 
radius of r = 0.01 by NUT lenses with N = 0.5 (solid curve) 
and A'' = 0.3 (dots), and by a Schwarzschild lens (circles) 



the stellar disk. When A^ = 0.5, the light from the 
source cannot be detected when u < 0.486. We es- 
timated the integrated fluxes (22) numerically with a 
relative uncertainty of 10~^. 

Figure 5 displays the V light curves for a star with 
r = 0.01, barely different from a point source, with time 
measured in units of the Einstein cone angular radius 
crossing time, Ie, for uq = 0.6 and to = 0. The curves 
are symmetrical about the maximum flux, and we show 
only their positive branches. The curve for N = 0.5 
(solid) is considerably different for the Schwarzschild- 
lens curve, plotted as circles. In agreement with (19), 
the observed differences between the lenses quickly de- 
crease with decreasing A''. The light curve for N ~ 0.3, 
plotted as dots, already displays little difference from 
the curve for a Schwarzschild lens, with the difference 
becoming virtually undetectable on the scale of the fig- 
ure when N = 0.1. No chromatic effects are observable 
for any N values in the studied range, < iV < 0.5. 
Deviations of the star's color indices from {B — V)o and 
{V — R)o arc within the expected calculation uncertain- 
ties. 

Figure 6 displays the results of our computations for 
a star with r ~ 0.10, comparable to the Einstein cone 
angular radius, for uq = 0.6 and to = 0. The results for 
lensing of the star by a NUT lens with A^ = 0.5 is shown 
by the solid curve, while the dots show the light curve 
for a point source lensed by the same lens. The curve 
with circles is the light curve for lensing of the same star 
by a Schwarzschild lens, and virtually coincides with 
that for the point source on the scale of the figure. We 
can see that, in contrast to the Schwarzschild lens, the 
influence of the star's finite size is fairly appreciable in 



Fig. 6 Light curves for microlensing a point source (dots) 
and a star with angular radius r = 0.10 (solid curve) by 
a NUT lens with A'^ = 0.5. Also shown is the light curve 
for microlensing of the same star by a Schwarzschild lens 
(circles) 



the case of the NUT lens. Thus, the difference between 
the Ivit) curves for the two lens types becomes more 
significant as the star's angular size increases. 

Figure 7 shows variations of the B — V and V — R 
color indices for the previous lensing case (?■ = 0.10, 
uo = 0.6, and to = 0). The solid curves show the 
rcsuhs for the NUT lens {N = 0.5) and the dots the 
results for a Schwarzschild lens. Despite the smallness 
of the chromatic effect, it is significant for the NUT lens, 
whereas the star's color indices during lensing by the 
Schwarzschild lens are virtually undistinguishablc from 
the initial values, {B — V)o and {V — R)o. Obviously, 
the chromatic effect can be made stronger by a special 
choice of narrow observed spectral intervals. However, 
patrol stellar microlensing observations use broadband 
systems. 

Another observable effect related to the predominant 
amplification of light from the disk limb can be man- 
ifest during stellar microlensing: linear polarization of 
the light can arise. The larger amplification coefficient 
of a NUT lens for fixed u in the source visibility zone 
enhances this effect compared to a Schwarzschild lens. 
On the other hand, it causes a shear of the source im- 
ages, rotating them in the plane of the sky; this should 
give rise to a rotation of the polarization plane, depolar- 
izing the combined light flux. The simple approach to 
estimating this effect used in section 3 is not applicable 
here, and this problem must be resolved by analyzing 
the polarization properties of stellar images computed 
using ray tracing. 
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Fig. 7 Variations of the B — V and V~R color indices dur- 
ing microlensing of a star with angular radius r = 0.10 by a 
NUT lens with A*' = 0.5 (solid curve) and by a Schwarzschild 
lens (dots) 



5 Conclusion 

We assume that small-scale clusters of weakly interact- 
ing massive particles can, in the course of their evo- 
lution, acquire rotational momentum and spheroidal 
shape. Even a small degree of oblateness similar to that 
of the critical surface for a Roche model can result in the 
appearance of caustics in the plane of a lensed source. 
The multiple images that form cannot be resolved us- 
ing current telescopes, and the cluster is manifest as a 
non- compact, spheroidal gravitational microlens. The 
crossing of caustics in the case of relative motion of the 
observer, particle cluster, and lensed star can produce a 
large variety of flux curves, including some observed in 
actual microlensing events that have been interpreted 
as manifestations of binary gravitational lenses. We 
must therefore bear in mind the possibility that at least 
some such events could in reality be associated with 
WIMP clusters. 

The unambiguous identification of such an exotic ob- 
ject as a wormhole from observations of gravitational 
microlensing of stars is a rather complicated problem. 
If a source with a small angular size does not cross the 
circular caustic, it is impossible to discriminate between 
gravitational lenses with positive and negative mass. 
Even for stars with relatively large angular diameters, 
deviations from the light curve for a Schwarzschild lens 
are of the order of the observational errors. Distinc- 



tive photometric features of microlensing by a worm- 
hole are a sharp decrease of the flux considerably be- 
low the initial level during the caustic crossing and a 
symmetric light curve. However, unavoidable gaps in 
the observations may prevent the detection of the sec- 
ond caustic crossing, and the shape of the light curve 
in the vicinity of the maximum is approximately the 
same as for microlensing by a caustic of an ordinary 
binary gravitational lens. When the angular size of the 
star is comparable to 0.10 of the Einstein cone radius 
for sufficiently large impact parameters, the observed 
light curves should resemble those for microlensing of 
a compact object surroundd by an extended gaseous 
envelope. This case can be easily distinguished from a 
negative-mass lens through an analysis of the chromatic 
and polarization effects arising during the lensing. 

Our results demonstrate also that, during gravita- 
tional microlensing of stars, under certain conditions, 
the influence of the star's finite size and the brightness 
distribution across the stellar disk can facilitate the de- 
tection of lensing objects whose local space-time is de- 
scribed with the NUT metric. Such effects are virtually 
absent if the stellar angular radius is much smaller than 
the Einstein cone angular radius. If the size of the star 
becomes order 0.10 of the Einstein radius, lensing by a 
NUT lens exhibits considerable differences from lensing 
by a Schwarzschild lens, in terms of both flux amplifi- 
cation and color-index variations. 
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